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Topological insulators are new quantum states with he- 
lical gapless edge or surface states inside the bulk band 
gap. These topological surface states are robust against 
the weak time-reversal invariant perturbations, such as 
lattice distortions and non-magnetic impurities. Recently 
a variety of topological insulators have been predicted by 
theories, and observed by experiments. First-principles 
calculations have been widely used to predict topologi- 
cal insulators with great success. 

In this review, we summarize the current progress in 
this field from the perspective of first-principles calcu- 
lations. First of all, the basic concepts of topological in- 
sulators and the frequently-used techniques within first- 
principles calculations are briefly introduced. Secondly, 
we summarize general methodologies to search for new 
topological insulators. In the last part, based on the band 
inversion picture first introduced in the context of HgTe, 
we classify topological insulators into three types with 
s-p, p-p and d-f, and discuss some representative exam- 
ples for each type. 




Surface states of topological insulator Bi2Se3 consists of a sin- 
gle Dirac cone, as obtained from first-principles calculations. 
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1 Introduction In two-dimensional electron systems 
at low temperature and the strong magnetic field, the Hall 
conductance a xy takes quantized values O, which proved 
to have a fundamental topological meaning! 2]. cr xy can 6 
be expressed as an integral of the first Chern number 
over the magnetic Brillouin zone. Recently quantum spin 
Hall(QSH) state was predicted and observed in CdTe/HgTe 
quantum wellOIH. In this system time-reversal symme- 
try(TRS) is present, and spin-orbit coupling(SOC) effect 
plays the role of Lorentz force in QH effect. The con- 
cept of QSH can be generalized to three-dimensional(3D) 



topological insulators with TRSEO. The electromagnetic 
response of a topological insulator is described by the topo- 
logical term of S e = (0/27r)(a/27r) / d 3 xdtE • B with 
7r where E and B are the external electromagnetic 
fields 0. This indicates the physically measurable and 
topologically non-trivial response, which opens the door 
for experiments and potential applications of topological 
insulators. 

Both 2D (QSH) and 3D topological insulators have in- 
teresting physical properties f7l[8l[9l[TOl . In this review we 
focus on 3D topological insulators with TRS. In this field, 
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an important task is to systematically search for all topo- 
logical insulators. In this process, first-principles calcula- 
tions played a crucial role. Up to now, most of topological 
insulators were predicted first by first-principles calcula- 
tions, and observed subsequently by experiments. 

2 Theories and Methods 

2.1 First-principles methods Density functional 
theory (DFT) is a formally exact theory based on the two 
Hohenberg-Kohn theorems (HK) |[TT1 . but the functional of 
the exchange and correlation interaction is unknown in 
Kohn-Sham(KS) equation lfT2l . In order to do numerical 
calculations, the local-density approximation(LDA) fT2t 
and Generalized Gradient Approximation(GGA) H 1 311 1 4l 
are usually used to approximate the exchange and corre- 
lation interaction in KS equation. Based on recent expe- 
riences, LDA and GGA work quite well for the study of 
topological insulators, because most topological insulators 
found to-date are weakly correlated electronic systems. 

As we know, the conventional LDA and GGA first- 
principles calculations tend to underestimate the band 
gap l 1511161 . However the band gap is directly related to the 
possibility of the band inversion which is the key topolog- 
ical property 0. For example, sometimes LDA and GGA 
predict a negative band gap, whereas the band gap is pos- 
itive in reality |[T7l . This can cause serious problem to pre- 
dict topological insulators. So it is necessary to improve the 
calculations of the energy gap. The most effective method 
to calculate the band bap is GW approximation (18). Sim- 
ply saying, GW approximation considers the Hartree- 
Fock self-energy interaction with the screening effect. 
Though the GW method has been used to study topological 
insulators (19), this method is very expensive. Besides the 
GW method, the modified Becke- Johnson exchange poten- 
tial together with LDA(MBJLDA), proposed by Tran and 
Blaha in 2009E1, costs as much as LDA and GGA, but 
it allows the band gap with the similar accuracy to GW's. 
MBJLDA potential can also recover LDA for the elec- 
tronic system with a constant charge density, and mimic 
the behavior of orbital-dependent potentials as well. 

LDA+UEB, LDA+DMFTd and LDA+GutzwillerO 
are employed to study to strongly correlated electronic sys- 
tems(d and f electrons), because LDA often fails for these 
systems. In strongly correlated electronic systems, the 
electrons are strongly localized, and have more features 
of atomic orbitals. This case requires proper treatment of 
atomic configurations and orbital dependence. Both LDA 
and GGA don't include the orbital-dependent feature, this 
is why they fail to describe strongly correlated electronic 
systems. Based on this understanding, all of LDA+U, 
LDA+DMFT and LDA+Gutzwiller include the orbital- 
dependent feature in different ways. For example, the 
on-site interaction is treated in a static Hartree mean-field 
manner in LDA+U method which is simplest and cheapest 
method. It is often used for strongly correlated systems, but 
it does not work well with intermediately correlated metal- 



lic systems. The self-energy of the LDA+DMFT method is 
obtained in a self-consistent way. Up to now LDA+DMFT 
is most accurate and reliable method, but its computational 
costs are high. LDA+Gutzwiller based on Gutzwiller vari- 
ational approach is recently developed. This method works 
well for intermediately correlated electronic systems, and 
it is cheaper than LDA+DMFT. Though it is still an open 
question how well these methods work on strongly corre- 
lated systems, it is true that these methods could reproduce 
some results of experiments, and that they can help to un- 
derstand some novel results in strong correlated electronic 
systems. 

2.2 Spin-orbit coupling Generally the SOC de- 
scribes the interaction of a particle's spin with its orbital 
motion. For example, in one atom, the interaction between 
one electron's spin and the magnetic field produced by its 
orbit around the nucleus can cause shifts in the electron's 
atomic energy levels, which is the typical SOC effect. In 
the non-relativistic limit, SOC Hamiltonian from relativis- 
tic Dirac equation is written as (24), 

i/soc = "4^ <7 ' PX(Wo) (1) 

where Ti is Planck's constant, mo is the mass of a free elec- 
tron, c is the velocity of light and a represents the Pauli 
spin matrices. H soc couples the potential Vq and the mo- 
mentum operator p together. 

In the case of the single atomic system Vq is spherically 
symmetric, H soc can be simplified, 

H soc = XL - a (2) 

where A is the strength of SOC interaction. L represents 
the angular moment. But in solid systems, Vo is the pe- 
riodic potential which can be very complex, so it is dif- 
ficult to exactly calculate SOC effect. For this case, SOC 
effect is usually calculated with a second- variational pro- 
cedure. SOC interaction is the key to the band topology, so 
all first-principles calculations to study topological insula- 
tors should be carried out with SOC. 

2.3 The criterion of topological insulators There 
are four Z2 invariants^; ^^2^3) for three dimensional 
topological insulators, first proposed by Fu, Kane and 
Mele0. When z/ = 1, states are strong topological insu- 
lators which have topologically protected gapless surface 
states consisting of odd Dirac cones. These surface state 
are robust against time-reversal-invariant(TRI) weak dis- 
orders. If uq = and at least one of ^1,2,3 is non zero, the 
corresponded states are weak topological insulators which 
have surface states with even Dirac cones on special sur- 
faces. We can simply consider weak topological insulators 
to be stacked by layered two-dimensional QSH states. In 
the presence of disorders, the surface states of weak topo- 
logical insulators can be destroyed. When all ^0,1,2,3 are 
zero, states are conventional insulators. 
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2.3.1 With the inversion symmetry The calculation 
of Z2 invariants is very simple for the compounds with 
the inversion symmetry. The formula of Z2 can be just ex- 
pressed with the parity values at the eight time-reversal- 
invariant moments(TRIMs) l25l , 

8 

(-ir = n* (3) 

1=1 

and 

(-If* = II *i=(nl,n2,n3) (4) 

rife=l;n J - ? £fc=0,l 

where 

N 

fii=l[ 6m(^) (5) 
m=l 

A 7 " is half of the number of occupied bands, and (2m(^i) 
is the parity eigenvalue of the 2 rath occupied energy band 
at TRIM K i= ( nin2n3 ) = 7j(nibi + n 2 b 2 + n 3 b 3 ) where 
bi,2,3 represent primitive reciprocal-lattice vectors. 

2.3.2 Without the inversion symmetry For the 
compounds without the inversion symmetry, Z2 invariants 
also can be calculated by the numerical method proposed 
by Fukui et al based on the tight-binding method |26) . 
Firstly, Z2 formula of QSH state can be expressed with the 
Berry connection and the Berry curvature, shown by Fu 
and Kane, 

Z2 — 77— [ (f A(k) - [ F(k)} mod 2 (6) 
27r Job- Jb- 

with 

A{k) = iZ n (u n (k)\V k u n (k)) andF(/c) = V k x A(k) 

(7) 

where B and dB indicate half of two-dimensional(2D) 
tori and its boundary, respectively. In order to do numer- 
ical calculations, equation © can directly be rewritten 
to its lattice version. Secondly, for 3D case, we can de- 
fine six 2D tori as Zo(k x ,k y ,0), Zi(k x ,k y ,ir), Yo(k x ,0,k z ), 
Yi(k X9 7r,k z ), Xo(0,k y ,k z ) and Xi(ir,k y ,k z ). We can calcu- 
late the Z2 based on equation © for each of these six tori, 
as zo, z±, yo , yi, xq and x\. The four Z2 invariants of topo- 
logical insulators are obtained by isq = xqx^, v\ = x n , 
V2 = 2/tt and v<$ = z n . Xiao et al. first successfully used 
these formulas to evaluate the Z2 invariants of half-Heusler 
compounds by first-principles calculations l27l . 

2.3.3 Adiabatic argument Sometimes it is not nec- 
essary to directly calculate Z2 for the compounds without 
inversion symmetry. One can start from an according com- 
pound with the inversion symmetry, and then adiabatically 
change this compound to that without inversion symme- 
try. If the energy gap does not close in an adiabatic pro- 
cess, the topological property will not change. For exam- 
ple, the space group of a-Sn is Fd3m(No. 227) and the 
inversion symmetry is held in this structure. We can eas- 
ily know a-Sn is topologically nontrivial from the parity 



calculations 1 25 ]. Based on the adiabatic argument, one can 
conclude HgTe is topologically nontrivial. 

2.3.4 Surface states Gapless surface states of topo- 
logical insulators must include the odd number of Dirac 
cones on one surface, and these surface states are robust 
against TRI weak disorders. So the calculation of surface 
states is another useful method to judge the band topol- 
ogy. The simplest way to calculate surface states is based 
on the free-standing structure. It is true that this is very 
powerful method to calculate surface states, but only for 
the compounds with the inversion symmetry and layered 
structure, such as, Bi, Sb, Bi2Se3 and so on. For exam- 
ple, if the compounds do not have the inversion symmetry, 
the polarization field might cause serious artificial effect, 
especially for the compounds with a small band gap. In ad- 
dition, if the compounds are not layered structure, the dan- 
gling bonds on the surface might cause a number of com- 
plex topologically trivial chemical surface states which can 
mix with topologically non-trivial ones. The topological 
surface states originate from the topological property of the 
bulk electronic structure. Though the details of these sur- 
face states can be modified by the special dangling bonds 
and the reconstruction of the electronic structure on the sur- 
face, we address that the topological feature don't change, 
such as, the odd number of Dirac cones. The calculation of 
the free-standing model also costs a lot, because the vac- 
uum layer should be thick enough in order to avoid the 
hybridization between the up and down surfaces. 

Besides the free-standing model, maximally local- 
ized Wannier function(MLWF) methods I28ll29l can be 
used to calculate the surface states 13011311 . Essentially the 
MLWF method is a tight-binding method, but the differ- 
ence from the conventional tight-binding method is that 
MLWF method can exactly reproduce the band structure 
of first-principles calculations. But it is not easy to obtain 
MLWFs, because the transformation from Bloch func- 
tions to Wannier functions is not unique due to the phase 
ambiguity of the Bloch functions used in first-principles 
calculations. Marzari and Vanderbilt reported an effective 
method to obtain MLWF by minimizing the spread func- 
tion ^Z n ((^ 2 ) — (r) 2 ) |28l . In order to calculate surface 
states, first we carry out the first-principles calculations 
for 3D bulk structure and then transform Bloch functions 
to MLWFs. At the same time the hopping parameters 
H mn (H) = (nO\H\mTV) between Wannier functions are 
obtained. At the next step, we use these hopping parame- 
ters to construct the hopping parameters of the correspond- 
ing semi-infinite structure, and then iterative method can 
be used to solve the surface Green's function, 

G^^e + irj) (8) 

Where n denotes the unit cell along the surface norm, 
and a is the Wannier orbital in the unit cell. The ML- 
WFs method can predict surface states well for layered 
compounds. For example, the calculated surface states of 
Bi2Se3 with MLWFs method agree well with the ones of 
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Figure 1 (a) and (b) are the band structure of HgTe by 
LDA and MBJLDA methods, respectively. i~6,7,8 represent 
the symmetry of energy levels at F point. The solid red 
circles indicate the projection of the s orbital of Hg. The 
LDA band structure shows the i~g — 17 — i~6 band sequence 
which is not correct, but MBJLDA can calculate the correct 
band sequence as Jg — i~6 — i~V. 



Angle Resolved Photoelectron Spectroscope(ARPES)ll30l 
[32l . Usually we don't expect to predict the exact dispersion 
of surface states, because this method does not include all 
complex situations on the surface. On the other hand, the 
surface states obtained from the MLWF method originate 
from the topological property of the bulk electronic struc- 
ture, so this is an ideal method to judge whether one com- 
pound is topologically non trivial or not. 

3 3-dimensional topological insulators After the 
initial discovery of the 2D topological insulator HgTe(3j 
|4), a number of 3D topological insulators are found with 
the great effort of theorists and experimentalists l7ll9l[T0l . 
In the following, we classify the topological insulators by 
the type of the band inversion, because the band inversion 
has a clear and general physical picture for most topolog- 
ical insulators. Up to now, there are three basic types of 
band inversions (s-p, p-p, p-f) in topological insulators dis- 
covered so far. In the following discussions, we will take 
some representative compounds as examples for each type 
of topological insulators. 

3.1 s-p type The most important s-p topological in- 
sulator is HgTe[3,4] which has the zinc-blende structure 
with the space group F43ra(No.216). Before HgTe was 
found to be a topologically non-trivial compound, it had 
been widely studies experimentally and theoretically l34l 
[35l[36l . Unlike other zinc-blende compounds, HgTe is a 



semiconductor with the symmetry-protected zero energy 
band gap . The Hg has occupied shallow 5d levels which 
tends to be delocalized, so Hg has a large effective positive 
charge in its core. The Hg's s level, which forms i~6 state 
in cubic symmetry, is pulled down below the Te's p levels 
which split into Jg and i~7, by this effective positive charge 
of Hg's core. Finally the energy level sequence at F point 
shows the i~g — i~6 — F7 order, which we call the s-p-type 
band inversion. If we define the energy gap AE, 

AE = E r6 - E r8 (9) 

where the Ep 6 and Ep 8 are the energy levels for Fq an d i~8 
at the r point. HgTe has a negative AE because the s-p- 
type band inversion, so it is well known as a negative gap 
semiconductor. 

The normal LDA and GGA can predict the band in- 
version between 7~6 and r$, but the exact band sequence of 
i~8 — i~6 — i~7 cannot be obtained (36) . The LDA band struc- 
ture with SOC shows the r^ — Tj — i~6 sequence, shown in 
filfila). As we addressed in the above, MBJLDA method 
can correct the error of LDA band structure. The band 
structure with the MBJLDA method is shown in figQIb), 
which perfectly shows the correct r$ — Tq — F7 sequence. 

Bernevig, Hughes and Zhang first identified the band 
inversion in HgTe to be the key ingredient of its topo- 
logically non-trivial behavior Q. Its topological invariant 
can also be obtained by an adiabatic argument ll25l . As we 
know, if we replace Hg and Te by the same atom in the 
zinc-blende structure, the crystal structure will change to 
the diamond structure with the inversion symmetry. Luck- 
ily, in nature grey tin has the diamond structure with space 
group ¥d3m, and it also is a semiconductor with a nega- 
tive energy gap AE due to the s level below the p level. 
Because grey tin holds the inversion symmetry, its parity 
values at all TRIMs can be easily calculated. It is worth to 
note that though grey tin is a zero band gap semiconductor, 
we still can define the topological property for all of its oc- 
cupied bands. Based on the formulas proposed by Fu and 
Kane, its Z2 invariants are calculated to be (1;000) which 
indicate topologically non trivial. Here the key is that the s 
and p at r point have opposite parity values. The occupied 
s state forms F^ , whereas p states form F^f and i~g~ . Tak- 
ing grey tin as the starting point, we can make a thought 
experiment to adiabatically change grey tin to HgTe. In 
this process, the negative gzp(AE) is never closed, which 
means the grey tin and HgTe have the same topological 
property. So HgTe prove to be topologically non trivial 
with Z2 invariant(l;000). Besides this adiabatic argument, 
HgTe's Z2 invariants can also be directly calculated by the 
numerical method addressed above. 

Similar to HgTe, there are a big family compounds 
known as half-Heusler materials(XYZ) l37l which in- 
clude more than 250 semiconductors and semimetals. 
Half-Heusler compounds have the face-centered cubic(fcc) 
structure sharing the same space group with HgTe. Y and Z 
form zinc-blende structure which is stuffed by X. Usually 
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Figure 2 (a) shows the crystal structure of chalcopyrite compounds(ABC2).(b) shows the energy gap AE for various 
chalc< ' irted. 
The L gical 
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Figure 3 The band structure of Bi2Se3 without (a) and with (b) SOC. The blue dashed represents Fermi level, (c), the 
evolution of the band sequence at F point starting from atomic levels. The three stage (I), (II) and (III) represent turning 
on chemical bonding, crystal field and SOC effects step by step. From Ref ll30l 
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Figure 4 (a), the calculated surface states of Bi2Se3 with MLWFs tight-binding method for a semi-infinite structure with 
the surface normal(l 1 1). The more red means more states, and the blue region indicates band gap without any states. The 
clear surface states with the linear dispersion at r can be seen in the band gap.(b), the ARPES result for Bi2Se3 along the 
T-M dirction. From Ref l30ll32l 



X and Y are transition metal or rare earth elements, and Z is 
a heavy element. The band structure of half-Heusler com- 
pounds at r point near the Fermi level is almost the same 
with that of HgTe case, s state forms r$, and p states split 
into i~7 and i~g. Some of half-Heusler compounds have the 
band sequence r$ — Tg — F-j and some others have the in- 
verted band sequence i~g — T§ — i~7. The interesting thing 
is that half-Heusler family were independently reported 
almost at the same time by three theory groups|38,39,27|. 
Besides the topological property, half-Heusler compounds 
are a class of multifunctional materials 140114 II . such as, 
superconductivity and magnetism, due to transition met- 
als and rare earth elementals. So half-Heusler compounds 
might be the best platform to study the Majorana fermion 
in topological superconductors (42|, dynamical axion field 
in topological anti-ferromagnetic phase O, and quantum 
anomalous Hall effect(QAH) in topological ferromagnetic 
phase l44l . Recently, some experiments, such as, ARPES 
and transports, already have been reported for half-Heusler 
compounds I45ll46ll47l . 

Generally due to the cubic symmetry, many topologi- 
cally non trivial compounds(HgTe and half-Heusler com- 
pounds) are zero gap semiconductors with Fermi level 
through /g level at r point, and a uniaxial strain is usu- 
ally needed to break the cubic symmetry in order to open a 
finite energy gap|48|. Feng et al reported that chalcopyrite 
structure can naturally break the cubic symmetry [[33]]. The 
chalcopyrite structure(ABC2) is the body-centered tetrag- 
onal structure with the space group I42d(No. 122), which 
could be regarded as a superlattice of two cubic zinc- 
blende unit cell AC and BC, seen in fi^2a). In essence, the 



unit cell of chalcopyrite are the double unit cell of HgTe 
with naturally breaking the cubic symmetry, and we ex- 
pect that these two class compounds might share the same 
topological property. Feng et al found that it is true that 
some materials with chalcopyrite structure are topological 
insulators, shown in fi^b). 

Besides the compounds talked about above, there are 
a lot of other s-p-type topological insulators, such as, 
/3-Ag 2 Te|49), KHgSb family El, NaaBiEIl CsPbCl 3 
family (52) and so on. 

3.2 p-p type Due to the simple surface states consist- 
ing of a single Dirac cone, Bi2Se3, Bi2Te3 and Sb2Te3 
compounds ll^Ql[32l[55l[56ll571l58H quickly became topolog- 
ical insulators extensively studied worldwide. Especially 
Bi2Se3 has a big energy gap of 0.3eV which is much larger 
than the energy scale of room temperature. These com- 
pounds share the layered structure with a five-atom layer, 
called the quintuple layer(QL), as the unit cell with the 
space group R3m(No. 166). Two equivalent Se atoms, two 
equivalent Bi atoms and a third Se atom are in each QL. 
The coupling is the chemical bonding between neighbor 
atomic layers whin one QL, but the van der Waals type, 
which is much weaker, between two QLs. It is worth to 
note that the inversion symmetry is held in the crystal struc- 
ture. 

In the following, we briefly introduce the basic elec- 
tronic structure of this family compounds by taking Bi2 Se3 
as an example. First of all, the band structure without SOC 
shows Bi2Se3 to be a narrow band gap insulator. Both the 
bottom of conduction band and the top of valence band are 
at r point, seen in fi^Ha). After SOC is turned on, the bot- 
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Figure 5 Schematic for the comparison of the surface 
states of (a) first-principles calculations! 31], (b) tight- 
binding calculations! 53] and (c) ARPES experiment (54|. 
From RefOH 



torn of conduction band is pulled down below the top of 
valence band, and an interaction gap opens at the crossing 
of valence band conduction bands, seen in fi^b). Based 
on the parity calculations, Z2 invariants of Bi2Se3 are cal- 
culated to be (1;000) which mean topologically non triv- 
ial. The key to be topological insulators is that the bottom 
of conduction band and the top of valence band have op- 
posite parity values. The schematic of the band sequence 
at r point clearly tells the band evolution starting from 
atomic levels with three stages, shown in fi^Oc). Because 
the s levels are much lower than p levels, we just start 
from the atomic p levels of Bi(6s 2 6p 3 ) and Se(4s 2 4p 4 ). 
At the stage(I), the bonding and anti-bounding effect be- 
tween Bi and Se atoms are considered. All the atomic or- 
bitals are recombined into P0~ yz , Pl^ yz and P2^ y ^ z 
where '0' represents the third Se, and '1', '2' represent Bi 
and the other two Se, respectively. '±' represents the par- 
ity values. At the stage (II), after the crystal field is turned 
on, p xyz levels will split into p xy and p z . The levels of 
P1+ and P2~ are nearest to Fermi level. At the stage (III), 
SOC effect is further introduced. P1+ becomes two degen- 
eracy levels(Pl+^), and P2~ becomes two degeneracy 

levels(P2~^) due to the time-reversal symmetry. SOC ef- 
fect tries to pull PI +^ down and push P2~^ up. Finally 
if SOC is strong enough, the p-p-type band inversion will 
happen between P^-ttl anc * ^ztl' 

Due to the layered structure with inversion symme- 
try, both the free-standing model and the tight-binding 
model based on MLWFs can be used to calculate surface 
states. Fi^JJa) shows the clear surface states of Bi2Se3 



with a single Dirac cone at P calculated by the MLWFs 
tight-binding model. Almost at the same time of Zhang 
et al.'s theory prediction l30lh Hasan group reported the 
the topologically non-trivial surface states of Bi2Se3 by the 
ARPES experiment! 32], shown in fi^Jb). Comparing the 
theory and experiment results, we have to agree that first- 
principles calculations can successfully predict topological 
insulators, including the details of surface states. Recently 
a lot of experimental studies of topological insulators are 
focusing on these compounds, because these compounds 
are easily to be grown by all kinds of experiments. 

The topological insulator Bii^Sb^O. 07 < x < 0.22) 
alloy also belongs to the p-p type 1251 . Bulk Bi and Sb 
share a rhombohedral R3m structure which holds the in- 
version symmetry, and they both are semimetals with some 
tiny Fermi pockets around the TRIM L and T points, but 
there is a direct gap at every k point through the whole 
Brillouin zone(BZ). So we can define an imaginary Fermi 
surface in the direct gap. Based on the parity calculations, 
we confirm that Bi is topologically trivial with Z2 (0;000), 
and that Sb is topologically non trivial with Z2 (1 ; 1 1 1). 
The key difference of Bi and Sb is the band sequence 
of the conduction and valence bands at three L points. 
For example, the conduction band of Bi is P s , and the 
valence band is L a where 'a/s' indicate the — /+ par- 
ity. Differently, these two bands switched with each other 
in Sb. After carefully comparing the band structure be- 
tween Bi and Sb, Fu and Kane predicted that the insu- 
late phase of Bii^Sb^ (0.07 < x < 0.22) alloy must 
be a topological insulator. Subsequently, Hasan group ob- 
served the topologically non-trivial property of Bii-^Sb^ 
by the ARPES experiment (54) . But the details of surface 
states don't agree with ones of tight-binding (53) and first- 
principles calculations (3D . The schematic of the differ- 
ence among these results are shown in fi^5l We can see 
that the ARPES result indicates three surface states J2i 2 3' 
but two surface state 2 are only found by tight-binding 
and first-principles calculations. Zhang etal. argued that 
the extra surface state might come from the imperfect 
surface, but this still is an open question up to now. 

Following Bi2Se3 family, a number of other p-p Bi- 
based topological insulators are predicted by theories and 
observed by experiments, such as, TlBiSe2 family (59), 
SnBi2Te4 and SnBi4Te7 family (60) , and so on. 

3.3 d-f type There is no clear evidence for the limit 
of the energy gap size for topological insulators. How 
could we find new topological insulators with bigger en- 
ergy gap(>0.3eV)? One possible way to enhance the SOC 
energy gap is to consider the cooperation of the SOC in- 
teraction and other effects, such as, the electron-electron 
correlation. In this idea, topological Kondo insulators were 
proposed, and SmB6 as an example was predicted to be 
a topological Kondo insulator (62l . Though due to 4f or- 
bitals SmBeis a strong correlated system, it only has a 
tiny energy gap. Recently Zhang et al. predicted AmN 
and PuTe family compounds are d and f topological in- 
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Figure 6 The band structure for AmN compound with 
U=0eV (a) and U=2.5eV (b) from the LDA+U method. 
The thickness of the band is corresponded to the projected 
weight of the d character of Am. In r — X direction, one 
band withd character clearly comes down to cross with the 
valence bands which mainly have the f character, (a) rep- 
resents semi metal without a band gap, but (b) represents a 
finite band gap. From Ref ll6Tl 



sulators with strong interaction ll6T1l . All AmN and PuTe 
family compounds have rock-salt crystal structure with the 
space group Fm3m(No. 225), and the inversion symme- 
try is also held in this structure. All these compounds have 
been well studied by theories and experiments, known as 
mixed valence materials. Here we take AmN as an example 
to understand the band structure. The configuration of ac- 
tinide Am is 5f 7 7s 2 6d°. The SOC interaction is strong than 
Hund's rule, so the f orbitals split into high energy J = 7/2 
and low energy J = 5/2 states. Approximately, in AmN, 
Am forms Am 3 + with the configuration 5^78° 6d°, the 
states of J = 5/2 should be fully occupied, and J = 7/2 
states are unoccupied. But due to the derealization of 5f in 
Am, 5f states partly hybridize with 6d states with neighbor 
Am atoms. 

In the fee crystal field, d orbitals first split into t<i g and 
e g states, and t<i g level goes down to cross 5f below Fermi 
level along F — X direction, shown in fi^6] The band in- 
version happens at three X points. If only LDA calcula- 
tions are used, the full energy gap cannot open through the 
whole BZ. After the electron correlation is introduced with 
LDA+U method, a band gap can open up with the proper 
correlation parameter U. We have to address that the elec- 
tron correlation U is found to enhance the SOC in these 
compounds. Because there are three TRIM X point in BZ, 
Z2 invariants of AmN must be topologically non trivial. 



Further more, our conclusion suggests all the mix- valence 
compounds with rock-salt structure must be topologically 
non-trivial. Especially, there transport experiments showed 
PuTe (63) has a big energy gap around 0.2eV, and this gap 
can be enhanced to 0.4eV with the pressure. Many of these 
f compounds host all kinds of magnetic phases, so they 
might open the opportunity to study QAH effect and dy- 
namic Axion field. 

4 Summary and outlook In this review, we first 
introduced widely-used techniques within first-principles 
calculations including LDA and GGA, GW and MB JLDA, 
LDA+U, LDA+DMFT and LDA+Gutzwiller methods, be- 
cause they play a crucial role on the field of topological 
insulators. Then the basic concepts of topological insula- 
tors and some useful methods to confirm the topological 
property are summarized. We classify topological insu- 
lators found to-date into three types as s-p, p-p and d-f 
based on the clear band inversion picture. For each type of 
topological insulators, we take several typical compounds 
as examples with talking about the electronic structure and 
the topological property. 

Though many topological insulators have been discov- 
ered, it is still important to find more with desired proper- 
ties. First of all, a big band gap is important for the applica- 
tion of surface states of topological insulators. Up to now 
the biggest band gap is around 0.3eV in Bi2Se3 compound. 
Secondly, the transport experiments to detect surface states 
are still very challenging [64, 65, 66 ]. One reason is that the 
quality of samples is not good enough with a low mobil- 
ity. Another reason is that Dirac cone always coexists with 
some bulk carriers. In order to overcome this barrier, on 
the one hand, experimentalists are trying to improve the 
quality of samples. On the other hand, it is important to 
find other new topological insulators with functional prop- 
erties. In addition, it is interesting to study the cooperation 
of the topological property with other phases, such as, su- 
perconductivity, magnetism and so on. We hope that this 
review can provide some guidance in the search. 
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